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SYSTEM AND METHOD FOR DETERMINING PROPAGATION 
CHARACTERISTICS OF PHOTONIC STRUCTURES 

TECHNICAL FIELD OF THE INVENTION 

[0001] The technique disclosed herein is directed, in 
general, to photonic structures (including waveguides) and, 
more specifically, to a system and method for determining 
propagation characteristics of such photonic structures. 

BACKGROUND OF THE INVENTION 

[0002] There is currently a great deal of interest in 
the light propagation characteristics of microstructured, 
or "holey," optical fiber waveguides with novel cross- 
sections consisting of holes surrounded by glass. The 
holes may be empty or filled with a material chosen to 
influence the propagation. The ability to vary the 
transverse geometry due to advances in fabrication 
technology, combined with the large index contrasts 
possible with such structures, give multiple new degrees of 
freedom that potentially enable designs radically different 
than were previously possible with standard fibers. 

[0003] Numerous interesting phenomena have already been 
observed in such waveguides. Among them are (i) guiding by 
the interference-based photonic band gap effect in fibers 
with an air core and a (truncated) transverse periodic 
lattice of air holes; (ii) variability of chromatic 
dispersion with microstructure ; and resulting (iii) 
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nonlinear effects in newly accessible spectral ranges due 
to microstructure- induced shifting zero dispersion point in 
glass core fibers. 

[0004] Efficient and accurate mathematical modeling of 
light propagation characteristics of microstructured fiber 
is necessary for their design and analysis. A feature of 
most such structures is that they are inherently leaky due 
to the existence of paths leading from the core to the 
cladding that avoid the holes and pass only through the 
background glass. Such structures support no true guided 
modes, however, they may have leaky modes characterized by 
complex-valued propagation constants or effective indexes 
(scattering resonances) ; the leakage rate is given by the 
imaginary part . 

[0005] Physically, this leakage is due to a combination 
of tunneling through the holes and propagating through the 
glass surrounding them. The ability to calculate such 
rates is clearly of fundamental importance. Other 
quantities of interest include the real parts of the 
effective indexes, which determine the response of the 
structure to longitudinal variations such as gratings, and 
the dispersion relations of the various leaky modes, which 
may be quite unusual compared to standard waveguides, due 
to the presence of large index contrasts and interference 
effects . 

[0006] Many numerical studies of microstructure fibers 
have been undertaken based on direct numerical calculation 
of static Maxwell's equations in an effort to determine the 



system's modes. Some of these numerical studies are also 
able to capture the attenuation rates. A variety of 
methods have been used, among them: (1) the multiple 
expansion, which works well for structures with circular 
holes, (2) more general expansions in local bases and 
Fourier decompositions, which are applicable to more 
general geometries and (3) scalar and vector beam 
propagation, which are applicable to general geometries but 
have limitations in computing very small attenuation rates 
and have proven problematic in some more complex 
geometries . 

[0007] All of these techniques have the characteristic 
that the computational difficulty of the calculations 
increases with the complexity of the structures. 
Accordingly, this creates a potentially intractable issue 
with today's proposed fiber structures. 

[0008] Accordingly, new systems and methods are needed 
in the art for mathematically modeling complex photonic 
structures that are more accurate, flexible, 
computationally tractable and free of the constraints 
prior-art methods. Such new systems and methods will allow 
the propagation characteristics of the photonic structures 
to be analyzed in detail, leading to the discovery of new 
photonic structures having highly advantageous 
characteristics . 
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SUMMARY OF THE INVENTION 



[0009] To address the above-discussed deficiencies of 
the prior art, the present invention provides a novel 
technique for studying microstructures and systems and 
methods designed around the technique. Systems for, and 
method of, determining propagation characteristics of a 
photonic structure having a transverse N-fold symmetry. In 
one embodiment, a system includes: (1) a numerical analyzer 
that employs a leading order systematic homogenization 
expansion having multiple scales to develop an angularly 
averaged indexed profile for the photonic structure and (2) 
a principal corrector, associated with the numerical 
analyzer, that employs details of the photonic structure 
and the homogenization expansion to obtain effective 
refractive indices of modes of the photonic structure. 
[0010] Further, in one embodiment of the present 
invention, a scalar approximation may be used to obtain 
accurate values for the real part of the effective indexes 
of the leaky modes for fine microstructures. For example, 
a multiple scale perturbation theory of a large class of 
"suf f iciently oscillatory" structures, suitable for the 
analytical study and efficient numerical computation of 
such quantities as leakage rates, group velocities and 
dispersion, that becomes more accurate as the transverse 
structure becomes more oscillatory may be used. 
[0011] If the transverse structure is invariant under 
rotation by 2nlN (e.gr., in the case of a single ring of N 
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holes uniformly distributed in an annulus) the leaky modes 
can be expanded, as well as the effective indexes 
(scattering resonances) in powers of UN . The leading order 
in the expansion turns out to be the averaged or 
homogenized index profile referred to above. The technique 
disclosed herein, thus, focuses primarily on the derivation 
and computation of the leading order corrections (order 
AT 2 ) , due to microstructure, or both the real and imaginary 
parts of the effective indexes of the leaky modes of such 
averaged profiles. Since the leading order behavior is 
given by a homogenization (in angle) effective medium, it 
may be referred to as a homogenization expansion. 
[0012] In an alternative embodiment, results on first- 
order corrections to homogenized eigenvalues of periodic 
composite media may be obtained in a different context. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0013] For a more complete understanding of the 
technique disclosed herein, reference is now made to the 
following descriptions taken in conjunction with the 
accompanying drawings, in which: 

[0014] FIGURE 1 illustrates a cross-section of a 
microstructured waveguide; 

[0015] FIGURE 2A illustrates a ring of N = 30 air holes 
approximated by a simple layered potential with 11 layers, 
wherein the ratio Rhoie/Rring of the hole radius to the ring 
radius is 0.7^/N; 

[0016] FIGURE 2B illustrates a graph depicting the 
averaged potential V av with R r ing = 12um; 

[0017] FIGURE 3A illustrates a ring of "air wedges" 
supported by N = 2 0 webs of glass, wherein the ratio Ri n /Rout 
of the inner to outer radii of the annulus is 10/11 and the 
air fill fraction within the annulus is / = 0.8; 

[0018] FIGURE 3B illustrates a graph depicting the 
averaged potential V av with Ri n = 10um; 

[0019] FIGURE 4A illustrates an 18-hole subset of a 
hexagonal lattice (N = 6) with inter-hole spacing A = 2 . 3um 
and hole radius Rhoie = 0.46 pm, approximated by a simple 
layered potential with L = 51 layers; 

[0020] FIGURE 4B illustrates the averaged potential V av 
in units of (pm) " 2 verus radius in units of pm; 
[0021] FIGURE 5 illustrates a block diagram of one 
embodiment of a system for determining propagation 
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characteristics of a photonic structure having a transverse 
N-fold symmetry constructed according to the principles of 
the present invention; 

[0022] FIGURE 6 illustrates a flow diagram of one 
embodiment of a method of determining propagation 
characteristics of a photonic structure having a transverse 
N-fold symmetry carried out according to the principles of 
the present invention; 

[0023] FIGURE 7 illustrates the effective indexes and 
attenuation coefficients of a set of the least-lossy 
resonances with /<e{0, 1, 2, 3} of the structure depicted in 
FIGURE 2, scaled to three different ring radii; 
[0024] FIGURE 8 illustrates the effective indexes and 
attenuation coefficients of a set of the least lossy 
resonances with /e{0, 1, 2, 3} of the structure depicted in 
FIGURE 3, for two choices of fill fraction x and three 
different radii; 

[0025] FIGURE 9 illustrates an attenuation of the lowest 
order (fundamental, or LP 0 i) resonance for a structure of 
the type shown in FIGURE 3, with R in = 1 jam, R out = 2 jam for 
the fill fractions / = 0.8, 0.9, and 1, with N = 3 and N = 

6 holes; 

[0026] FIGURE 10A illustrates results of the weak 
residual calculation |/^ , (tf)-f£r* ,l, (tf)| and (0)|| ; 

[0027] FIGURE 10B illustrates results of the weak 
residual calculation N 2 * (6) - F^ k ' x) (0)) for different N; and 

[0028] FIGURE 11 illustrates an example of a simple 
layered approximation to a structure containing six 
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circular air holes. 
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DETAILED DESCRIPTION 



[0029] A few of the consequences of the homogenization 
expansion and numerical implementation that the technique 
disclosed herein wishes to highlight include: 

[0030] (1) Computational algorithm: The analytic theory 
leads to a natural efficient computational algorithm for 
the modes and spectral characteristics, e.g., leakage rates 
and dispersion. Multiscale analysis enables one to 

eliminate the "stiff" aspects of the computation due to the 
rapidly varying structure. Therefore, one can expect 
significant reduction in computational effort. This is 
especially important for increasingly microstructured 
media, for which the approximation improves. This, 
however, is also especially important for simulation-based 
optimization of light-guiding characteristics. 

[0031] (2) Arbitrary geometry and index contrasts: 

Although the technique disclosed herein requires N-fold 
symmetry of the structure, the individual microf eatures may 
have arbitrary geometry and large index contrasts. In one 
implementation of the technique disclosed herein, an 
arbitrary microfeature is approximated by a simple layered 
potential, defined below, for which many of the 
calculations required can be done explicitly. 

[0032] (3) Very good agreement with full numerical 

simulation: Results of numerical simulations are presented 
below based on a theory for several waveguides with 
transverse microstructure . The technique disclosed herein 
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further focuses on the complex effective indexes as they 

vary with the geometry of the microstructure . The 

technique disclosed herein, in one particular embodiment, 

is particularly interested in the imaginary parts of the 

effective indexes, which correspond to the leakage rates. 

The technique disclosed herein may then present comparisons 

of predicted loss rates, with results obtained by solving 

the Helmholtz equation directly by Fourier decomposition 

algorithm and the vector Maxwell equation by multiple 

methods. One expects a homogenization theory to be valid 

when the wavelength of light is long compared with the 

individual microf eatures . The technique disclosed herein 

finds very good agreement with the theory using full 
simulations even in regimes where X fs ld , the ratio of 

wavelength to length scale of microf eatures , is as small as 

3/2 . Although the expansion is derived for N being large, 

with the N = 00 limit being the leading order term, examples 

set forth below show agreement in cases where N = 3 and N = 6. 

The technique disclosed herein also notices the expected 

departure of the approximate methods from the results of 
full simulations for sufficiently small k f Jd . 

[0033] (4) Sensitivity of leakage rates to 
microstructure: The technique disclosed herein uses a 
theory to compute the first two nontrivial terms of the 
effective indexes (scattering resonances) of the leaky 
modes. The first term corresponds to an average theory and 
the second term is a correction due to the microstructures . 
The imaginary parts of the effective indexes, corresponding 
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to the leakage rates, are very sensitive to the 
introduction of microstructure and their accurate 
approximation requires both terms. In contrast, as noted 
above, the real parts of the effective indexes may be 
relatively insensitive to microstructure, and are well 
captured by the leading order term. 

[0034] (5) Corrected fields: Corrections, due to 

microstructure, of the mode fields predicted by the 
averaged structure are compactly supported in space if the 
microstructure perturbation is compactly supported in 
space. This is typically not possible for the true 
solution of the mathematical model. An additional 

observation indicating a limitation of the homogenization 
expansion is seen for a family of structures consisting of 
two rings of holes, parameterized by the relative phase 
angle of the arrangement in one ring to another. Since the 
variation in refractive index in one ring occurs on a set 
disjoint from that for the second ring, the homogenization 
expansion predicts complex effective indexes that are 
independent of the phase angle. This suggests that in the 
regime where certain interference effects are important, 
the expansion technique described herein may have limited 
use. Nevertheless, a two- term approximate expansion of 
modes and effective indexes proves valid for sufficiently 
large N . This rigorous theory gives field corrections 
which, at any finite N, are not compactly supported. 

[0035] Mathematically, the leaky modes and effective 
indexes are obtainable from the solution of the eigenvalue 
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problem for the Schrodinger equation 



(1) 

y/ outgoing as r = \x\ -> oo 

[0036] This is formulated as a scattering resonance 
problem below. Here, A ± denotes the two-dimensional 

Laplace operator in the transverse plane. (When the 
context is clear, the subscript 1 will be omitted.) Since 
the outgoing radiation condition at infinity is not self- 
adjoint, one expects the eigenvalues, E , to be complex. 
The attenuation or leakage rates are given by the imaginary 
part of 



P = y \k 1 n 2 g -E (2) 

where n g denotes a background refractive index. Of 
importance in optics is the effective index 



n~ s — 

eff k 



(3) 



where k = 2nlX fs , and X fs is the free space wavelength. The 

effective index is a complex quantity for leaky waveguides. 
The real and imaginary parts of p , 9ty? and 3/? , are 

displayed in the FIGURES in plots of 
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Power Input _ n 

y = i01og 3£ 

Power Output 



versus 



Attenuation rates, /, are given in units of dBlcm , while <K„, # 
is dimensionless . 

[0037] The technique disclosed herein, in one aspect, 
also relates the certain portions of the invention to a 
rigorous perturbation theory of a general class of 
scattering resonance problems with rapidly varying and high 
contrast potentials. In particular, this perturbation 
theory implies the validity of the homogenization expansion 
and error estimates for the two-term truncation uses in 
certain simulations . 

[0038] Maxwell's Equations in a Waveguide 

[0039] Consider a waveguide with refractive index 
profile n(x). If the medium is non-magnetic and there are 
no sources of free charge or current, then Maxwell's 
equations imply that the transverse components of the 
electric field obey 

where ( 4 ) 

v x e L =-V 1 (e x -V 1 lnn 2 ). 



[0040] In one embodiment, the technique disclosed herein 
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works in the scalar approximation, which entails neglecting 
the term v x e L in equation (4) that arises due to the vector 

nature of the fields. For modes that are localized away 
from the microstructure, which includes most of the 
examples given below, the vector corrections are expected 
to be small because the field e L is small at the interfaces, 

which are the only places at which it contributes. 
However, this will not be true in other examples of 
interest, and one skilled in the art would expect that the 
methods can be extended to the full vector case. In any 
case, vector effects will not be discussed any further. 
[0041] In the scalar approximation the components of e ± 

satisfy independent scalar Helmholtz equations. For 
purposes of the technique disclosed herein, <p denotes 

either of these transverse components. Then, 



where A ± denotes the Laplace operator in the transverse 
variables, x ± =(x u x 2 ). Introducing the notation: 



the equation for the cp and propagation constant can be 
viewed as a Schrodinger equation with potential V(* L ) and 

energy parameter E : 



& L +k 2 n z )<p = p 2 <p 



(5) 



V = k 2 (n 2 g -n 2 ),E = k 2 n 2 g -p 2 



(6) 



L(p = (-A ± +V)(p = Eg). 



(7) 
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[0042] Corresponding to the physical problem described 
below, it can be assumed that for some n>0, 

K(x)>0,for all x and 

(8) 

K(x) = 0,for |*| >r.. 

[0043] The potential, V , therefore does not support 
bound states (guided modes) and only has scattering states 
(radiation modes) , along with scattering resonances (leaky 
modes) . 

[0044] The resonance problem for micro structure 
waveguides 

[0045] Now, consider a cylindrical waveguide composed of 
two materials w g" and tt h ,; , with corresponding refractive 
indexes n g and n h . One can specify a transverse cross- 
section of the waveguide in the 2 -dimensional plane, with 
coordinates x = (x,y) = (r,0) . In FIGURE 1, shown is a schematic 

picture of a glass waveguide 100 with a single ring of 12 

air holes 110 in a solid glass cylinder 120. Four regions 
are illustrated: A x the region |x|</? lVl , a disc of glass with 

refractive index n g ; A 2 , the region R in < \x\ < R out , an annular 

region with 12 equally spaced air-holes; A 3 ; the region 
R out < \x\ < R clad , consisting of a solid glass n = n g , and finally; 

A 4 , the region \x\ >R clad . consisting of air, n = n h . 

[0046] The technique disclosed herein, in part, is 

interested in how well light is confined to the core (Ai) 
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region in such microstructure fiber. This structure has 

actual guided, in addition to leaky, modes, due to the 

lower refractive index of region A 4 . However, the 

propagation of light from Ai to A 3 may be considered as 

loss. Therefore, the appropriate mathematical problem for 

these purposes is the mode equation (7) posed on the 
spatial domain {|jc| < /? to }v^{/^ <|jc|</? OH/ }u{/? OItf <|x|} and subject to an 

outward going radiation condition as |x|-»oo. in practice, 

the glass optical fiber is usually coated with a high index 
and high loss polymer, so any light that does leak into A 3 
is rapidly attenuated, further justifying the use of this 
boundary condition. 

[0047] A more general mathematical formulation of the 
problem is now presented. 

[0048] The scattering resonance problem: Given a 
potential (index profile) v = V(r 9 0) of compact support, the 

technique disclosed herein seeks E and nonzero <p for which 
L<p=(-A + V)<p = E<p (9) 

subject to the outward going radiation condition, as r->ao . 
The pair (<p,E) may be called a scattering resonance pair. A 

scattering resonance pair is sometimes also called a 
quasiresonance . 

[0049] In this setting, the radiation condition can be 
simply described. Because V(r) is identically zero for 
r>r m , equation (7) becomes: 
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-(A x +£)p = 0, r*r.. 



(10) 



[0050] Solutions of equation (10) can be expanded in a 
series 




(11) 



where ffj l) and h\ 2) are first and second Hankel functions of 

order / . Therefore the requirement that <p be outgoing is 

equivalent to the vanishing of all coefficients cf ] , in 
equation (11) : 

Outward going radiation o c| 2) = 0 for all /. (12) 

[0051] Since equation (11) is a non-self adjoint 
boundary condition at infinity, E can be expected to be 
complex. However, there are constraints on the location of 
E in the complex plane. For the <p to be oscillatory at 
infinity, <K£>0, and for q> to be outgoing, W£>0. 

Furthermore, 3£<0. Also, if 3£>0 then 3V£>0. From the 
series expansion (equation (11)) with c, (2) =0, it can be 
concluded that <p exponentially decays as r-»oo and therefore 
is an eigenstate of the self-adjoint operator L with a non- 
real eigenvalue. This yields a contradiction. 
[0052] Thus, if E is an energy for which a solution of 
the resonance problem exists, then WE>Q and 3£<0. 
Correspondingly, by (2.3), ^>0 and 3>3>0. Since 
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E(x,a>) = e(x ± )e ifix * , the complex resonances E correspond to the 
radiative decay or leakage rates of the structure. 



be derived to the resonance problem for microstructures, 
giving a homogenized (averaged) theory at leading order, 
plus systematically computable corrections. The analysis 
is carried out for a class of potentials V with the 
dependence : 



where V is 2d. periodic in □ and □. Thus equation (13) 

allows both a slow and a fast angular modulation of the 

index profile. In the special case where V does not depend 
on □ and V = V(r,N0), the potential corresponds to an index 

profile with an N-fold symmetry. Assume V>0 and that for 
some r, s V(r) = 0 forr > a; . The technique disclosed herein will 

show the following: 

[0055] Homogenization Expansion: The resonance mode 
problem (equations (9) through (12)) has, for large N, 
solutions with the formal expansion: 



[0053] 



Multiple scales and homogenization expansion 



[0054] 



A multiple scale expansion of solutions will now 



V = V(r,0,N0) = V(r,0,@) 



(13) 



(p(r,e,N) = e {N) (r,0,e), e = N0 



given by 



(14) 




E (A,) = 
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Thus, corrections due to microstructure begin at o(jV~ 2 ) . 
(<D o (r,0),£ o ) is a non-trivial solution of the resonance 

problem: 



L av Q> 0 =E 0 ® 0 (15) 



where O 0 is subject to the outward going radiation 

condition (equation (12)), and the averaged operator is 
given by: 

L av ^± 2 ]Lde^-A^V av (r,e) (16) 



* O = o (p) + o (/ ° , where <D (P) is the mean zero solution (easily 

2 2 2 2 •* 

derived from the, Fourier series of K(.,.,0) of: 



p-^* ( 2 p) =[V(T,0,0)-VJr,0)]%(r 9 0) (17) 
and (o ( 2 A) ,£ 2 ) solves: 

[0056] Finally, E 2 is uniquely determined by the 
condition that (4.6) has a solution that satisfies an 
outgoing radiation condition at r = oo; see equation (12). 
[0057] This yields an approximate solution of the scalar 
wave equation and an approximation (neglecting vector 



effects) of the transverse electric field, E = (E 1> E lM ) = (E i> E 2 ,E 3 ) , 
of Maxwell's equations: 



E 



(x-p) ~ e^^[* 0 (\xJ[Aa>) + ±O 2 ^AWa) + O 



Here, q = 0 or 1, co = ck , and f3 = y Jk 2 n 2 g -E =tt/? + /3/?, with3jff >0 . 
E ±9 (jt x ,jt 3 ) decays with increasing jc 3 and is therefore a "leaky 

mode." These modes are not square integrable. 
[0058] The results of numerical simulations performed 
using a two-term truncation of the homogenization expansion 
will be presented and compared with published computations 
of solutions to the full problem. Further the question of 
convergence of the expansion is discussed. 
[0059] A. Multiple scale expansion 

[0060] Due to the rapidly varying coefficient in 
equation (7) , one skilled in the pertinent art would expect 
rapid variations in its solutions. Therefore, the 

technique disclosed herein explicitly introduces the fast 
angular variable: 

® = N0 (19) 

and views <p as a function of the three independent 
variables : r,0 and 0 : 

<p(r,0) = O(r 9 e,®). (20) 
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Here, the polar coordinates are x = (r,0) and note that 

r 



where A r denotes the radial part of the transverse 
Laplacian: 



A =3 2 +-3 . 

r v r r' 

r 



Thus, the operator d e is replaced by d 0 + Nd Q and equation (7) 
becomes : 

~(d 0 + Nd Q ) 2 + K(r,tf,0)j* = E®. (22) 
Equivalently 



2 A 

0 = £0. (23) 



Seeking an expansion of O and £ in the small parameter l/N 



0 N 1 N 2 N 3 N 4 5 

E (N) =E + 1 E + _L E +E (H) 

0 N 1 N 2 3 



(24) 



Here, O y . = O(r,0,©) . 

[0061] Substitution of equation (24) into equation (23) 
and equating like powers of l/N, a hierarchy of equations 
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arising at each order in l/N is obtained. The first five 
equations of this hierarchy are: 

o(Af 2 )4^(D 0 =0 (25) 



0(N):^ 9 <b l ^~d 9 d e <b 0 (26) 



0(1) : \d 2 Q <S> 2 = ~d § dj> x +(L- E 0 )O 0 (27) 



O(^ 1 ):^a^ 3 =-43A ( I>2+(L-Eo)^-E 1 0 0 (28) 



0{N->)±d^=-^d e d Q e^ (29) 



[0062] B. Construction of the multiple scale expansion 
[0063] As is typical in perturbation expansions, a 
hierarchy of inhomogeneous linear equations, which have a 
fixed linear operator to invert and varying inhomogeneous 
term, must be solved. Each member of the hierarchy is of 
the form: 



\d 2 Q F = G(r,d,Q). (30) 
r 



A necessary and sufficient condition for solvability in the 
space of 2n-periodic in 0 functions is: 

\G{r,e,p)dp = 0. (31) 
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Proceeding with a term-by-term construction of the 
perturbation expansion: 

o(n 2 ) Terms : 



-1-^*0*0 (32) 
r 



From this it follows that <D 0 is independent of the fast 
phase, 0 : 

* o M,©) = * o (r,0) (33) 
O(N) Terms: 



-A d e d e®o ~3e*i=0. (34) 
r r 



Using equation (33) in equation (34) 



(35) 



and therefore: 



0,(^5,0) = 4), (r,tf). (3 6) 



Using equation (36) , the following can be found: 
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0(1) Terms : 



-^o 2 =(l-e 0 )o 0 . (37) 



[0064] Equation (37) has a solution in the space of 2n- 
periodic in 0 functions if and only if the following 
solvability condition holds: 

(4 v -£ 0 )O 0 =0. (38) 

Here L av and V av denote, respectively, the average operator 
and potential with respect to the fast angular dependence. 
They are given by: 



(39) 

1 f 2<r 



Note that L-L av = V -V m . 



[0065] (<D o (r,0),£ o ) are fixed to be a scattering resonance 

pair of an outgoing solution and its associated complex 
eigenvalue . 

[0066] Since <D 0 satisfies equation (38) , the 
inhomogeneous equation (37) can be put in the simpler form: 

-^<t 2 =[v(r,e,e)-VJr,e)]<!> o (r,0). (4 0) 



The solution of equation (40) is expressed in the form 



(41) 



where O ( 2 p) (r,0,0) andQ™(r,0) are, respectively, a particular 

solution of mean zero and a homogeneous solution of (40) , 

which is to be determined at higher order in N" 1 . That 
d> 2 p) (r,0,Q) can be chosen to be 2 ;r -periodic in 0 follows from 

the mean zero property of K(r,0,0)-K av (r,0) . Also, a>^ } is mean 

zero in 0 because the constant term in its Fourier 
expansion is included in 0 ( 2 A) . To compute <£ ( 2 p) , V(r,0,®)-V av (r,0) 

is expanded into a Fourier series in 0 : 

V(r 9 O 9 e)-VJr 9 ff) = Zti(r 9 0) e » (42) 

Resulting in 

Q>[ p) (r,0,®) =d; 2 [F(r,0,0)-F (r,0)]r 2 O o (r,0) 

(43) 

0{N' 1 ) Terms: 

\d 2 ® = -4-5 a <D + (L-E )0> -E* 

(44) 

The solvability condition becomes 
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which is satisfied by taking £, =0 andO, =0 . 

Using the mean zero property of <J> ( 2 p) , one can see that 
equation (44) becomes 

<5 e O(r,<9,0) = -2d e O[ p) (r,O 9 0) (46) 



0(AT 2 ) Terms: 

Using that £, =0 andO^ =0 , as well as equation (46) to solve 
for 0 3 in terms of 0 ( 2 p) , it can be found that this order: 

\d 2 ® =4-3 2 <I> (p) +(L-E )<D -EO 

r 2 e 4 y.2 9 2 V (K 2 2 0 

= (L-E +-id 2 )0> (p) +( L -E )<D {y,) -E O />ir7 ^ 

V 0 r 2 fi y 2 V 0 7 2 2 0 (47) 

= (L -E +4"^)0 (p) +(L-L W A) 

v av 0 j.i 9' J v av' 2 

= (L -E#-EO +(L-L )<J> (p) . 

v av 0' 2 2 0 v av 7 2 



[0067] The first two terms on the right hand side of 
(4.35) have zero average in 0. Therefore, the solvability 
condition for D 4 becomes: 



(L -E )O w =EO -— f (L-L )<& {p) dp 

v av 0 / 2 2 0 2ft Jo av y 2 r 



= E 2 0> o --L \ 2 y(r,e,p)-VJr > e)]^[ p) (r,0,p)dp (48 ) 



where equation (43) is used and where <p~ p y denotes the mean- 
zero antiderivative . Thus, an outgoing solution of 
equation (48) is sought. This determines E 2 and therewith 

the expansion for the scattering resonance (effective 
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index) through order N' 1 . 

[0068] This truncated expansion will now be used 
numerically to obtain the complex effective indexes 
(scattering resonances) for various structures. The 
details of the implementation are set out below. In 
particular, an explicit expression for E 2 (see equations 

(76), (77) and (69)) for the class of AT -fold symmetric 
waveguides is derived. 

[0069] FIGURE 5 illustrates a block diagram of one 
embodiment of a system, generally designated 500, for 
determining propagation characteristics of a photonic 
structure having a transverse N-fold symmetry constructed 
according to the principles of the present invention. 
[0070] The system 500 includes a numerical analyzer 510. 
The numerical analyzer 510 employs a leading order 
systematic homogenization expansion having multiple scales 
to develop an angularly averaged indexed profile for the 
photonic structure in the previously described manner. The 
system 500 further includes a leading order corrector 520. 
The principal corrector 520. The principal corrector 520 
is associated with the numerical analyzer 510 and employs 
details of the photonic structure and the homogenization 
expansion to obtain effective refractive indices of modes 
of the photonic structure, again as previously described. 
[0071] FIGURE 6 illustrates a flow diagram of one 
embodiment of a method, generally designated 600, of 
determining propagation characteristics of a photonic 
structure having a transverse N-fold symmetry carried out 
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according to the principles of the present invention. 
[0072] The method 600 begins with a start step 610, when 
it is desired to analyze a photonic structure to determine 
one or more of its physical characteristics. The method 
600 proceeds to a step 620 in which a leading order 
systematic homogenization expansion having multiple scales 
is employed to develop an angularly averaged indexed 
profile for the photonic structure in the previously 
described manner. The method 600 then proceeds to a step 
630 in which details details of the photonic structure and 
the homogenization expansion are employed to obtain 
effective refractive indices of modes of the photonic 
structure, again as previously described. The method ends 
in an end step 64 0. 

[0073] Numerical Simulations for Selected Structures 
[0074] The technique disclosed herein now illustrates 
the theory set forth above with numerical calculations 
performed for several classes of structures. The details 
of the implementation are set out below. The three classes 
of structures considered are: 

[0075] Example 1: a ring of circular air holes 200 
(FIGURE 2 A) , 

[0076] Example 2: an annulus of air supported with glass 
"webs" (FIGURE 3A) and 

[0077] Example 3: a subset of a hexagonal lattice 
(FIGURE 4 A) . 

[0078] The technique disclosed herein will identify 
individual leaky modes or scattering resonances, of these 
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structures using the LPi m (linearly polarized) notation 

appropriate for solutions of the scalar wave equation. 

Usually this notation applies to guided modes, but no 

ambiguity will result, as none of the structures considered 

support -guided modes. The subscript 1<={0, 1, 2,...} refers 

to the angular dependence e iie of solutions to the averaged 

equation, while me{l, 2, 3 ... } indexes the collection of 

leaky modes with fixed 1. This collection is ordered by 
SRn^ , with m = 1 corresponding to the largest such value 

(when this rule is applied in guided modes, the usual 
meaning of m as the number of nodes in the radial 
wavefunction is recovered) . 

[0079] In the FIGURES, distinct plotting symbols are 

used to represent different values of 1; the various values 

of m for fixed 1 are not explicitly displayed but can be 
inferred from the monotonic relationship between m and Wn^ . 

[0080] Examples 1 and 2: The geometries of Examples 1 

and 2 were varied to produce a total of 3 0 different 

structures between them: The ring of holes was sealed to 

three different core sizes, while the annulus with webs was 

scaled to three different radii for each of two different 

air-fill fractions and three different periodicities 2^/N. 

Each of the structures was approximated by a simple layered 

potential. (The web-supported annulus depicted in FIGURE 

3A is a simple layered potential without approximation) . 

In all of the examples, the refractive index of glass was 

taken to be 1.45, and the holes were assumed empty. The 
calculations were all performed at X f5 =1.55//m . All of these 
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structures support only continuous spectrum, so the complex 
eigenvalues associated with resonances of the structures 
are reported on. 

[0081] The following sets of resonances for the rings of 

circles and webbed structures depicted in FIGURES 2A and 3A 

have been computed. Attention has been restricted to leaky 

mode solutions of the averaged structure with effective 
indexes Wn^ € [1 . 41 , 1 . 45] , attenuation coefficients less than 

about 2dB/mm / and angular index le{0, 1, 2, 3}. To express 
these quantities in a notation consistent with the 
technique disclosed herein, recall from equation (6) that 



£ = J*V-E (49) 

and that the effective index is given by 

« =k' l fi. (50) 

Also, the attenuation y (in dB/cm, if p has units of l//zm) 
of a leaky mode is given by 



y = 3/? 51) 

r ln(10) H 



[0082] The factor of 2 in equation (51) is due to the 
definition of / as the attenuation coefficient for 

intensity, rather than field amplitude. Once this set of 
leaky modes of the averaged potential is found, the leading 
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order corrections E 2 to the leading order energies E 0 can be 
computed. 

[0083] The results for the three-scaled rings of circles 
(FIGURE 2A) are presented in FIGURE 7. Also, those for the 
three scalings and two air-fill fractions, along with 
various values of the periodicity N , of the ring of "air 
wedges'' are presented in FIGURE 8. 

[0084] Resonances of a structure of the type shown in 

FIGURE 3 with R in = 1/xm, R out = 2/xm for the fill fractions f 

= 0.8, 0.9, and 1, with N = 3 and N = 6 holes were also 

computed. The calculations were performed for a range of 
free-space wavelengths X fs ranging from 1/im to 2/im. The 

structure with f = 1 is an idealized ring of air with no 
supporting structure. In FIGURE 9, the results of the 
theory to parallel calculations performed by a Fourier 
decomposition algorithm were compared. It can be observed 
that the two methods agree exactly for the f = 1 case, as 
they should. 

[0085] Additionally, in all cases the present method 
agrees with the Fourier calculation as A, fs ->2{im. Even for 

smaller X fs the agreement is quite good except for the f = 

0.8, N = 3 case, which is not surprising given the small 
value of N and the fact that the width of the ring is equal 
to the smallest free space wavelength considered. The 
plots displayed in FIGURE 9 are consistent with the 
expectation that approximation by the homogenization 
expansion improves (1) for fixed k fs and increasing N, as 

well as (2) for fixed N and increasing X fs . 
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[0086] The dashed lines in FIGURE 9 display the 
attenuations of the averaged structures, without the o(w~ 2 ) 

correction. The necessity of including those corrections 
is evident, though the effect in these structures is not as 
dramatic as in some of the other structures considered, 
including that of Example 3 below. 

[0087] Another result is the effective index of the 

first excited state above the fundamental (the LPn state) 
for the f = 0.9, N = 3 structure at X fs = 1.55/xm. In Table 

1, the results of the Fourier method with truncations of 

the homogenization expansion to one and two terms, for both 

the fundamental and first excited leaky modes are compared. 

It can be noted that the one-term (averaged theory) 

truncation predicts the real part of the effective index 
well, but again, the o(at 2 ) corrections are necessary to 

obtain good agreement of the attenuation rates. 



Mode 


0 tn Order Homogenization 


2 nd Order Homogenization 


Fourier 




Neff 


Atten 


n e ff 


Atten 


n eff 


Atten 


LPoi 


1. 3712 + 0. 0000387i 


1.36 


1. 3727 + 0. 0000644i 


2.27 


1.365+000071i 


2.48 




1. 2468+0. 000436i 


15.3 


1.2515+0. 000832i 


29.3 


1. 255 + 0. 00075i 
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Table 1 - Homogenization Expansion Versus Fourier 



[0088] Example 3: Finally, the structure depicted in 
FIGURE 4 is considered. This is an 18 -hole subset of a 
hexagonal lattice (N = 6) with inter-hole spacing A = 2.3/zm 
and hole radius Rhoie = . 46//m. Notice that the fundamental 
(LPqi) resonance has a leakage rate of 14 dB/cm, while that 
of the averaged structure was 0.92 dB/cm. By comparison, 
the solution to the full vector problem using a multipole 
method with outgoing radiation conditions results in a rate 
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of 16dB/cm. In this example, the effect of the N~ 2 E 2 
microstructure correction on the leakage rates is even more 
apparent than in Example 2. However, it cannot be said how 
much of the discrepancy is due to vector effects, and how 
much to the other approximations in this method. 
[0089] Convergence as N->co 

[0090] A consistency check on the implementation of the 
scheme described above for simple layered potentials will 
now be described. The first non- trivial corrections N 2 ® 2 
and N~ 2 E 2 are computed to the field and energy, 
respectively, in equation (24) and express the field and 
the approximation plus error terms: 



<D = 0 +-^-<X> +0 (N) 
0 N 2 3 

E = E +-i-E +E (N) . 

0 N 2 2 3 



(52) 



[0091] It is then considered by how much the mode 
equation (9) is not satisfied by the first two terms of the 
expansion equation (52) . Therefore, the residual is 
defined as: 

Fra ^_ A+K _ ( E o+ -^ E )]. (53) 



Substituting equation (52) into equation (9) , it is 
observed that the residual is given by 



F =E (N) 



1 



S(N) 



$ + — o +<x>' 



A + M E + — r-E 

0 N 2 2 



(N) 



(54) 
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[0092] F res is expected to be of order N" 1 , because 
and are formally of order AT 3 and the Laplacian in 

(54) , when acting on the fast dependence, gives a factor of 

N 2 . Indeed, in expanding F res in powers of 1/N, it is found 
that when 1 > 0, the leading term, , is the term 

F^>=-Lai*, . (55) 

res fly* 0 3 



This may be simplified using equations (46) and (43) to 
FM=-l(dO)( d -i(V-V )). 

res Jy|\ 0 0j\ 0 \ av ) ) 



[0093] Therefore, when 1 > 0 it is expected that as w->oo 
F ^F (l) ~of-landF -F<" 1) ~ of-U . (57) 

res res 1 N J res res 1 ^2 I 



when f = o(<X> o =a> o (r)J, d g ® o vanishes and the leading behavior of 
the residual is higher order in l/N. 

[0094] Since the potential V or its derivatives may have 
discontinuities, a problem is encountered when attempting 
to compute F res using equation (54) . For example, it is 
possible that 0 ( 2 p) is discontinuous as a function of r. 

Equation (54) is, therefore, interpreted as holding in the 

"weak sense" or sense of distributions, and both sides are 
integrated against a radial test function <j> R (r) of compact 

support. This allows one to move the radial part of the 
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Laplacian over to act on $ R and to define the "weak 
residual" 



Fi^>(*) = ]nfr^;-^ (58) 



and its leading behavior in 1/N 



(59) 



[0095] Equations (58) and (59) have been computed for 

the structure of Example 2 (FIGURE 3) with R in = 10/zm, R out = 

11/xm for the LPn leaky mode. Further, N is allowed to vary 

from 20 to 22,000. The results are presented in FIGURE 8. 
The solid line in FIGURE 8A represents \\^ k \0)- F^ kA) (0)\\ m as a 

function of N. The scales are logarithmic on both axes so 
the slope of -2 verifies that equation (57) holds over the 
entire range considered. The dashed line in FIGURE 10A 



represents 




It should be noted that the slope 



changes from -2 to -1 at around N = 1000, which again 
verifies equation (57) . 

[0096] FIGURE 10B shows the pointwise behavior of the 
weak residuals by plotting N 2 n(F}™ k \0)-F^ i] (0)) vs. 0 for N = 20 

and N = 54. The envelope of the residuals is invariant 
after scaling by N 2 , suggesting that the N~ 2 term dominates 
the residual after the loading term is subtracted, even for 
those moderate values of N. This behavior persists for all 
values of N considered, up to N = 22,000. 
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[0097] Implementation of Homogenization Expansion 

[0098] The procedure used to obtain the simulation 

results above will now be outlined. It is specialized to 

the case of a potential with no slow angular dependence, so 
V(r,0,®) = V(r > ®) . In this case the potential has an N-fold 

symmetry. 

[0099] This discussion pertains to N-fold symmetric 
microstructure, whose individual microf eatures have 
arbitrary geometry. The strategy is to then approximate a 
general microstructure potential V by a simple layered 
potential. For the case of simple layered potentials, 
analytic expressions can be obtained which facilitate the 
numerical computations . 
[0100] General Structures 

[0101] The angularly average potential is a radial 
function, V av = V av (r) . One can take 

O o (r,0) = /(r)^ 
/ being an integer. Then, f (r) satisfies 

(^ + 7 + U')-b.)/M-o (60) 

[0102] A solution of the resonance problem is a pair 
(f(r),E Q ) , such that /(r) is a nonsingular solutions of 

equation (60) satisfying an outgoing radiation condition at 
r = oo, equation (12). As noted above, E 0 is complex with 
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3£ 0 < 0 . 

[0103] In general, the existence of solutions to the 

resonance problem is a subtle technical problem. It will 

now be analyzed in detail for the particular class of 

simple layered potentials, which will be defined. The 
existence of a solution, (/(r),£ 0 ) will also be firmly set. 

Note that since V av (r) = 0 for r>r. , 

f(r) = kHU(JTr) i r>r (61) 

for some constant k, . 

[0104] The leading order correction, Ar 2 0 2 , is then 

calculated by the perturbation scheme. In the separable 
case, the equation for <X> ( 2 A) is of the form: 

( L „-*.X-(veWK (62) 

where Q(r) has compact support. Since O o (r,0) = /(r)e ,w , all of 
the Fourier modes of <D ( 2 A) (r,0) except the first vanish, and so 
<& i f ) (r,0) = U(r)e il9 is defined, and one can consider the equation 

■ (63) 

[0105] Recall that by equation (8) both Q and V av are 
zero for r>r. and one gets 



A+J-^HO-VW.'*'.- (64) 



38 



One can proceed with the solution of the inhomogeneous 
problem (equation (63)) for r>0 and then obtain an 
expression for E 2 by imposing the outward going radiation 
condition. 

[0106] It has been assumed that V is regular at the 
origin. One can let y^(r) denote the solution to the 

homogeneous equation associated with equation (63) that is 
regular at r = 0, and denote an independent solution, which 
together with % spans the solution set. Then J>,(r)ocr' and 

y 2 ( r ) ocr ~' if />0 > or j?i(r)oclnrif/ = 0 . For r>r, , an appropriate set of 

homogeneous solutions of equation (64) is in terms of 
Hankel functions, h\ j) {JlT 0 r}j = \,2 . Hf ] satisfies the outward 

going radiation condition at infinity and Ff} 2) the inward 
going radiation condition at infinity. 

[0107] A particular solution of equation (63) can be 

constructed by variation of parameters using these 
homogeneous solutions in the intervals R in =(0,r.) and R out =(r„cc). 
g in [F] is defined to denote a choice of particular solution 

of 



g 0U ,[F] to denote a choice of particular solution of 



-A +-=■-£„ 

r r 1 0 



gjF] = Hr)> r M reR o»Y 



(66) 
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Then, one gets 



inV U(r)=E 2gii [f](r) + gjQf](r) 



(67) 



in R m : U{r) = (^r) + ^ (^r) + ^ [/](r) 



(68) 



where the coefficients 4 and rj, are determined by imposing 
continuity of U and d r U at r = r, : 



2/y 



=v " (r '\vK[/]W-^[/K0)J 



kffllto 



(69) 



with £ and 77, given by equation (69) , equations (67) and 

(68) define a nonsingular solution of equation (63) . It 

remains to impose radiation condition at r = oo . Consider 
U(r) for r >r, . Let 



[0108] Then, equation (64) becomes: 



z 2 G"+zG , + (z 2 -/ 2 )C7=^-*z 2 // ( (,) (z), r(z)>r. (70) 
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By equations (68) and (61) 



G(z) = t /(r(z))^(£j//;')( Z ) + ,.(£j<)(z)3^[//; i >](z). (71) 



g out can be constructed by variation of parameters. A 
particular solution is 

^(i)jf(»)+^(*)i5fW 

where 



7j {p) \h](z) = + r , ' v A d$ 



(72) 
(73) 



Therefore, 



Using equation (74) in equation (71), one gets: 

G( Z ) = t/(r(z))^(£)i/f)(z) + J7 .(£)^(z) 

> (Of ^ ^ (i) (75) 

0 
0 

Recall that £(^2) an d *7.(E 2 ) are determined by the continuity 
conditions at r = r. , and are given by equation (69). The 
condition determining E 2 is that the incoming part of G(z) 
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be identically zero. Derivation of this condition on E 2 
requires the use of certain integral identifies involving 
Hankel functions. 

[0109] Note that both terms in equation (75) that are 
proportional to are outgoing, as can be seen by 

referring to their asymptotic forms. Now one gets 

f4^(<r)] 2 *r -£([^(0j-^(0^(fl) • 



[0110] Therefore, the terms in equation (75) that are 
proportional to Hj 2) can be written as: 



([^)(z)] J -^ i (z)^(^))^ ) (4 



Hf\ 2 ) 



ikn £ z 2 
4 £ 2 

0 



Again from the asymptotic form of H\ x) , the latter term is 
seen to be outgoing at infinity so the condition E 2 is: 

a(^)-^4(W(«.)T-«S(«.)«S(«.))-o (76) 



7. (£ 2 ) can then be written as 

n.(E t )^+E 2 ft. (77) 



where rj, 0l can be read off from the second component of 
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equation (69) . Equations (76) and (77) give a linear 
relation for E 2 that can be solved. 

[0111] B. Simple layered potentials: Leading order 
resonances 

[0112] An implementation of the above scheme is now 

outlined for a case of a simple layered potential, which is 
defined to be one with potential V(r,0,&) that is both 

independent of the slow angular variable 0 and is a simple 
function of the form 

rM)-^^MW e > K « (78) 

where L is the number of radial layers and Mi is the number 
of angular sectors in one period of the i th layer. 1^(£) 

denotes the indicator function of the set A, taking the 
value one for £ e A and zero for £<£A. By definition, 
r x = 0, 0, s 0 and 0^, = 2/r . The potential in equation (78) has, as a 
function of 0, M t jumps in layer Ri and, as a function of 
0NM i , jumps in layer Ri. An example of a simple layered 
approximation to a structure containing six circular air 
holes is presented in FIGURE 11. 

[0113] With the definition in equation (78) of V(r,e), V av (r) 
is constant on each region Ri = [r±, r i+1 ] ; equivalently, 

V (r) = V reR. (79) 

av v ' av,t i 

In the special case where the index profile, r(r,6) , takes on 
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two values rig and n h and 2nf i denotes the total angle of the 
annulus Ri occupied by "material h" , one gets by averaging V 
in equation (6) : 



v J r )= kl f,[ n r<l r * R ,- (80) 

[0114] In each internal Ri, the radial wavefunction f(r) 
can be expressed as a linear combination of Bessel and 
modified Bessel functions: 

/«-V.(*( r )) + Vi(*( r )) 

where the definitions of y lf y 2 and z depend on the sign of 
% K- K fl v) :z = ( ± K- K av )) 1/ V; fori<L,y r J i orI r y 2 =Y i orK; while in the 

outer region i = Lwe let y x =H?\ y2-H\ 2) . The condition that the 

solution be nonsingular at r = 0 implies that 



CHI) 



(82) 



or a multiple thereof. The values of (o;.,r.) in some interval 
i determines the values (cr ul ,T M )by continuity of f (r) and its 

derivative at the interface. This relationship may be 
expressed by a transfer matrix Ti(E 0 ), which gives rise to 
the relationship between the coefficients in the outermost 
layer and those in the innermost : 



u 
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where 



T h)- T W KKK) T ,W 



(84) 



[0115] 



The radiation condition is equivalent to the 



orthogonality of (0, 1) and (g l ,t l ) . Therefore, the resonance 
or scattering energies are determined by the following 
scalar transcendental equation: 



where T k j , denotes the (k, j) energy of the matrix T(E 0 ) . 

C. Simple layered potentials: Microstructure corrections 

to $> 0 and E 0 

The order N ~ 2 correction to <D 0 will now be constructed, 
which is N" 2 0> 0 and to E 0 , which is N ~ 2 E 2 for the special case 

of a simple layered potential given in equation (78) . 
[0116] Recall that <D 2 =0 ( 2 p) +0 ( 2 A) . Here ®[ p) is, in general, 

given by the Fourier series of equation (43) . When the 

potential is of the form given in equation (78) , 0 ( 2 p) can be 

explicitly calculated as described below. solves 

equation (48) and E 2 is chosen so that equation (43) has a 



(o i)t(e o ) = 0 



(85) 



or, equivalently 



T 2] (E o ) = 0 



(86) 
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solution satisfying the outgoing radiation condition at 
r = oo. We now focus on the determination of E 2 and <£ { 2 A) in the 
separable case. Recall that the equation (48), for 0 ( 2 A) , is 
of the form: 

( L .-*>?=( £ , + G( r )h (87) 

where Q(r) has compact support; its explicit form for 
simple structures is presented below. Since ® 0 (r,&) = f(r)e ilB , 

one can define O ( 2 A) (r,0) = U(r)e ae and consider the equation 

[0117] The solution of equation (88) in region R q can be 
written as 

U W = U l P) (*(')) + V, (»(r)) (89) 

where the definitions of y lt y 2/ and z are the same as in 
equation (81) , and U { q p) denotes a particular solution 

obtained by variation of parameters. As above, the 
coefficients (<^ +1 ,77 9+1 ) can be obtained from (£ 9 ,^) by imposing 

continuity of U(r) and its derivative at the interface. 
[0118] The details are omitted, but note that the 
transfer matrix formulation from equation (82) is modified 
in this case to reflect the presence of a homogeneous term: 
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(4 ^ 



*7 



= F (E ) + T 



(90) 



Note that the matrix T q does not depend on E 2 , while F q 
does . 

[0119] The solution in the innermost region R X/ Ui (r) , 
must be regular at the origin. Through a choice of 
particular solution, one may take the homogeneous part of 
the solution in this region to be identically zero. Then, 
iteration of equation (90) yields: 



= F +S F +S S F +..S S ..SF, 

L I 1-1 I L-l L-2 L i-l 2 1 



(91) 



The radiation condition in the outermost region R L remains 
to be imposed. Note that £?(r) = 0 in (7.28), and that by 

construction D 0 is outgoing and thus can be written as 



(92) 



[0120] Therefore to impose the radiation condition on 
the correction and thereby obtain E 2/ equation (70) with 
k. = <* L may be used. The calculations above can therefore be 

used to obtain a particular solution. This gives the 
equation for E 2 : 



(93) 
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where 



7] (E ) = tj +Er] (94) 

is read off from equation (91) . E 2 is then found by solving 
the linear relation given in equations (93) and (94) . 
[0121] Calculations for Simple layered potentials 
[0122] Expressions needed to compute the leading 
corrections to O 0 and E 0 due to microstructure will now be 

presented. Two quantities are required: 

0<; } (r, 0, 0) = r 2 0 (r, 6) 6£ [v (r, 0) -Vjr)] (95) 



and 



[0123] The mean-zero antiderivatives in equations (95) 
and (96) are explicitly calculable due to the form of the 
potential (78). For reJ^ , 

^^^-^Wl-cr^e+l^ce^+^je-ej) (97) 



where 



Sj~ir u (* M -* t ) (98) 
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and the constant Ci is given below. Integrating again, once 
more for re^., 



5- [V (r. 0) - K (r )] = C, + C ,0 - 1 ^ ( © 2 + 



(99) 



where 



(100) 



and the constants Ci and C 2 are given by 



c-«c + WLv - 

2 1 £ av,i 



(101) 



— vfrf© -0 U-S (© -0 f +-P (©. -©.A 

Finally, for re^,., the integral in equation (96) evaluates 
to 

if|a;[K(r.p)-K p (r)]|»*-C' + ic i r j|W + 

KvU'K.- 0 ,) 3 } 



[0124] Conclusion 
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[0125] From the above, it is apparent that the technique 
disclosed herein provides a systematic homogenization 
expansion for microstructured cylindrical waveguides with 
transverse N-fold symmetry. The method of derivation and 
multiple scale analysis facilitates the removal of fast 
scales due the rapid variations in the microstructure . 
Therefore, numerical implementation of the expansion does 
not encounter the intrinsic "stiffness" associated with 
problems having large separations of scales. The technique 
disclosed herein has further implemented the expansion for 
a number of sample structures and has computed the 
effective indexes of their leaky modes (scattering 
resonances) . Of great importance are the imaginary parts 
of the effective indexes, which correspond to the leakage 
rates and result from a combination of propagation and 
tunneling losses. In contrast to the real parts of the 
effective indexes, these leakage rates are very sensitive 
to the geometry of the microstructure in the case of low 
index (air) holes. 
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